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Multiplicities, fluctuations and correlations* 
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The recent results on hadron multiplicities in heavy and light quark fragmentation, multiplicity local fluctuations 
and multiparticle correlations submitted to the Conference are reviewed. 



1. Introduction 

One of the most important observables in par- 
ticle production processes is multiplicity, i.e. the 
number of particles (mostly, hadrons) produced 
in the collision [1]. The multiplicity dependence 
on the energy scale, species of particles, event 
flavour content are among the main predictions of 
the theory of strong interactions, quantum chro- 
modynamics (QCD) [1,2]. On the other hand, 
the multiplicity is used to select or to describe 
events, e.g. as a trigger for specific processes, as 
an input for kinematic variables' spectra. The 
distribution of multiplicity, its mean value and 
multiplicity fluctuations are the essential charac- 
teristics of the collision dynamics. However, the 
multiplicity distribution tells us just about the 
averaged, integrated numbers, while deeper in- 
formation comes from the moments of the distri- 
bution, which measure particle correlations, i.e. 
probe the interaction dynamics [3] . 

Here, I report on the multiplicity flavour 
dependence'^, on the analyses of local multiplic- 
ity fluctuations and multiparticle correlations'*. 
These studies provide us with details of strong 
interactions and allow us to estimate the level of 
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the applicability of QCD, based on the partonic 
picture, to the production of hadrons which are 
the experimentally observed objects. 

2. Definitions and notations [1,3] 

The multiplicity distribution, or the density 
Pm of multiplicity n of particles with kinematic 
variables pi,p2, ■ ■ ■ ,Pn is defined by the inclusive 
probability spectrum, 

1 dn{pi,p2, ■■■,Pn) 
Ncv dpidp2 ■ ■ ■ dpn 



Pn{Pl,P2, ■ ■ • ,Pn) 



where A'ov is the number of events. 

As it follows from this formula, the single parti- 
cle distribution pi{pi) gives an average multiplic- 
ity, (n) = J pi{pi)dpi, while integration of the 
g-particle densities leads to the unnormalised qth 
order factorial moments, 

fl=j Pl(Pl'P2 • ■ ■Pn)dpidp2 ■ ■ ■ dpn 

= (n(n-l)---(n-g+l)) = (nM). (1) 

The normalised factorial moments are then de- 
fined as Fq — fq/{n)'^. 

The g-particle densities give us a way to 
study particle correlations described by the q- 
particle correlation functions, (factorial) cumu- 
lants, Cq{pi, . . . ,pq). The cumulants vanish 
whenever one of their arguments is statistically 
independent, i.e. these functions measure gen- 
uine g-particle correlations. 

The cumulants are constructed from the mul- 
tiplicity densities, e.g. 

Ci{pi) = Pi(pi), 
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C2{Pl,P2) = P2(Pl,P2) - Pl(Pl)Pl(P2), 
C3(Pl,P2,P3) = P3{Pl,P2,P3)- 

^Pl{Pl)P2{P2,P3) +'2pi{pi)pi{p2)Pl{P3) ■ (2) 
(3) 

These functions, being properly normalised, are 
used to study muZtiparticle correlations in differ- 
ent kinematic variables. 

In studies of local fluctuations and correlations, 
i.e those in phase-space bins, one uses the bin- 
averaged factorial moments and cumulants. The 
phase space is divided into M bins of equal size, 
so that {rijn) is the number of particles in the mth 
bin. Then, the normalised bin-averaged factorial 
moment is defined via Eq. (1) averaged over bins. 

These factorial moments have been exten- 
sively used to extract the non- statistical (non- 
Poissonian) fluctuations in many types of colli- 
sions [3] . Such fluctuations lead to the power-law 
scaling of factorial moments as a function of M 
called the intermittency phenomenon. 

In analogy with Eq. (3), Eqs. (2) are averaged 
to the bin-averaged normalised cumulants, 

KJM) = — y ^P-^, (4) 

where fc^™' are the Mueller moments, the func- 
tions of unnormalised moments fq{m), Eq. (1), 

fci = (n™), = 
fcf^ = (^.[?J)-3(n[2l}(n„,) +2(n„,}^ (5) 

3. Hadronisation of heavy and light quarks 

In this Section, I consider the recent results on 
hadron multiplicities from fragmentation of heavy 
and light quarks [6-8]. The study of the quark 
content in multiparticle production provides one 
of the basic tests of QCD. The results from LEP 
are of a special interest since they cover a wide 
centre-of-mass energy region and can be directly 
compared with QCD which is mostly predictable 
at the asymptotic energies [2]. 
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Figure 1. The difference in the mean charged 
multiplicities, Shi, between heavy and light quark 
pairs as a function of centre-of-mass energy. The 
dashed-dotted line is the combined result from all 
measurements, see Ref. [8] for more details. 

In Refs. [6] and [8], a study of the fragmen- 
tation of heavy b-quark and light quarks (1 = 
u, d, s) is performed. The measurements of the 
difference in charge particle multiplicities, Shi = 
(n^h) ^ {''hl)^ bb and 11 events in e"'"e~ anni- 
hilation at the centre-of-mass energies above the 
peak are carried out. The findings are com- 
pared to the theoretical predictions of QCD and 
to a more phenomenological (the so-called naive) 
model (for review, see Ref. [2]). The QCD calcu- 
lations predict energy-independent behaviour of 
the multiplicity difference Shi, while in the naive 
model one expects the decrease with increasing 
energy. The latter is connected with the assump- 
tion that the hadron multiplicity accompanying 
the heavy hadrons in bb events is the same as the 
multiplicity in 11 events at the energy left to the 
system once the heavy quarks have fragmented. 

The difference between the heavy and light 
mean quark-pair multiplicities obtained by DEL- 
PHI at 206 GeV is Su = 4.50 ± 1.05(stat) ± 
0.52(syst) [8], while OPAL finds Su = 3.44 ± 
0.40(stat) ± 0.79(syst) in the energy range of 130 
- 206 GeV [8]. The difference in the values is 
connected with some differences in the data pro- 
cessing procedure. In the meantime, the results of 
the energy-dependence is found to be the same for 
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the both studies: the mean multiphcity difference 
is energy-independent and favours the QCD pre- 
dictions while it is inconsistent with the flavour- 
independent naive model (Fig. 1). Also in agree- 
ment with QCD is the ratio between light quark 
multiplicities, {rii) : (nj) ^ 1, where i,j ={u,d,s}, 
as obtained by OPAL at the peak energy [7]. 
To note is that the (riu) and (rid) are highly sta- 
tistically anti-correlated (—90%) which is due to 
their fractions in K-mesons. 

4. Scaling of local fluctuations in hadronic 
Z° decays 

During the last decades, the phase-space lo- 
cal multiplicity fluctuations are actively studied 
in many reactions, from leptonic to nuclear colli- 
sions, and the intermittency scaling of the fluctu- 
ations has been established [3]. 

All the studies show that the intermittency is 
more pronoimced in high dimensions, while in one 
dimension (e.g., in rapidity) the effect is diluted 
by projection onto one dimension. This leads to 
flattenning of the factorial moments. Such a be- 
haviour is well understood from a QCD parton 
shower which is a three-dimensional branching 
and naturally leads to the fractality. The best 
area to study the effect and its connection with 
QCD is given by e+e~ annihilation, where such 
investigations has been performed earlier [3,9,10] 
and new analysis is carried out by L3 [11]. 

The L3 analysis gives us further hints about the 
intermittency origin, which is currently far from 
understanding. The new study employs the fact 
that so far the fluctuations in many dimensions 
were studied with the same mimber M of bins 
in each dimension/direction. This leads to self- 
similar fractals, i.e. the fluctuations in any direc- 
tion are assumed to be the same, or isotropic. 
In case when the dynamics in different direc- 
tions are not equivalent, the fluctuations become 
anisotropic and this tells us about the scli- affinity 
[12]. A self-afflne behaviour has been observed 
in hadronic interactions [12], while a self-similar 
scaling has been qualitatively confirmed in e+e" 
annihilation [10]. In terms of factorial moments, 
Eq.(3), the 3-dimensional moments are expected 
to exhibit the intermittency property in e"'"e~ col- 
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Figure 2. The 1-dimensional 2nd-order factorial 
moments, Eq. (3), as a function of the number M 
of bins in the randomised and the thrust frames 
compared to the Monte Carlo predictions with 
and without Bose-Einstein correlations. The fit 
is given by the first f-la of Eqs. (7). The be- 
haviour in the qq frame is similar to that in the 
randomised frame. See text and [11] for details. 

lisions when phase-space is partitioned isotropi- 
cally, while this scaling in hadronic interactions 
only occurs for anisotropic partitioning. 

In order to quantitatively study the observed 
isotropic fluctuations, the method of the Hurst 
exponent is used [12]. The Hurst exponent, 

Hab = lnMa/\nMb, (6) 

is obtained from the fit of the one-dimensional 
second-order factorial moments, 

F^iMa) =Aa+Ba M--^" = (7) 

1 + 7a 

Here a, h are the directions of the (a, b) plane. For 
the isotropic dynamical fluctuations, Hab = 1, 
while Hab ^ 1 if the fluctuations are anisotropic. 

L3 uses rapidity y, azimuthal angle (p and 
transverse momentum pt - the variables often 
used in multiparticle studies - in the analysis. 
The variables are defined with respect to the 
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Figure 3. The all-charge and like-sign cumulants, 
Eq. (4), as a function of the number M of bins 
in rapidity and in three dimensions. The like- 
sign cumulants are compared to the Monte Carlo 
predictions with and without Bose-Einstein cor- 
relations. See text and Ref. [13] for more details. 

thrust axis, as well as the other two frames are 
considered: the randomised and the qq frames. In 
the former one, the (^-angles are randomly cho- 
sen being uniformly distributed in [0, 27r], while 
in the latter case one uses Monte Carlo to correct 
the original thrust axis to the qq axis. 

The 3-dimensional factorial moments are found 
to exhibit the intermittency scaling when phase- 
space is isotropically partitioned. This confirms 
earlier indications for isotropic fluctuations in 
e+e" annihilations [9,10]. No sensitive differences 
in the different frames as well as no disagreements 
between the data and Monte Carlos are observed. 

Fig. 2 shows the 2nd-order factorial moments 
in one dimension, from which the quantitative 
test of the isotropy of the fluctuations is made 
using Eqs. (7). The calculations of the Hurst ex- 
ponent give: Hab « 1, a, 6 = {y,ip,pt}, i.e. due 
to the definition of Eq. (6) one concludes about 
the equivalence of the directions and about the 



isotropy of the fluctuations. The fluctuations in 
If are nearly absent in the thrust frame. In all 
frames the data are well described by Jetset, while 
somewhat less well by Herwig. No sensitivity to 
Bose-Einstein correlations is seen. From its study, 
L3 also concludes that there is the dependence on 
the QCD dynamics which serves to decrease the 
fluctuations in the thrust frame. 

The 2-jet sub-samples are also analysed using 
the Durham jet algorithm. The fluctuations in 
the 2-jet events show the self-afflne (anisotropic) 
fluctuations. This looks like the hard gluon emis- 
sion leads to the isotropy. 

5. Multiparticle correlations 

The multiparticle correlation studies are per- 
formed for e"'"e~ annihilations at the Z° energy 
[13] and for pp annihilation at the incident mo- 
mentum of 22.4 GeV/c [14]. 

Genuine multiparticle correlations in e~'"e~ an- 
nihilations are studied by OPAL [13]. To extract 
the correlations, the method of the bin-averaged 
normalised cumulants is used, see Eqs. (4)-(5). 
In addition to the earlier analysis [10] of the cu- 
mulants of all charged particles ("all-charge cu- 
mulants"), OPAL studies the cumulants of same- 
charge particles ("like-sign cumulants"). The in- 
vestigation is performed in rapidity y, azimuthal 
angle f and transverse momentum pt, the same 
variables as in the above described L3 analysis, 
but the sphericity axis is used as a reference axis. 

Fig. 3 shows the all-charge and like-sign cumu- 
lants calculated in one (rapidity) and three di- 
mensions. Even in one dimension, positive gen- 
uine correlations of groups of two, three, four par- 
ticles are present. The cumulants exhibit a scal- 
ing behaviour, although in rapidity the saturation 
already appears at the moderate bin sizes. The 
like-sign cumulants increase faster and drive the 
all-charge ones at small bins while the unlike-sign 
cumulants are almost constant. This points to 
the likely influence of Bose-Einstein correlations. 

The comparison between the data and Monte 
Carlo like-sign cumulants shows that the model 
describes well the data when Bose-Einstein cor- 
relations are implemented. The same is true for 
the all-charge cumulants (not shown here). This 
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Figure 4. The Ochs-Wosiek plot in two and 
three dimensions for aU-charge and hke-sign cu- 
mulants. The dashed hne shows the Negative Bi- 
nomial prediction. The solid line is a fit, InKq = 
ttq + Tq \nK2. See text and [13] for more details. 

is not the case of the L3 results where no differ- 
ence is found for the models with and without 
Bose-Einstein correlations, see Fig. 2. To note 
is that the OPAL data is described well in one 
dimension as well as in three dimensions. 

From Fig. 3, one can see that Bose-Einstein 
correlations which are implemented in the model 
as the correlations of two identical particles and 
are pair-wise adjusted, well also describe the cu- 
mulants oi q > 2. This suggests to consider the 
interdependence of the 2nd- and higher-order cu- 
mulants. Fig. 4. The fit akin to the Ochs-Wosiek 
one [3] gives the same parameter values for the all- 
charge and like-sign cumulants, while disagrees 
with the Negative Binomial predictions. 

The multiparticle fluctuations are also studied 
in [14] using the Serpukhov fixed target pp an- 
nihilation data at 22.4 GeV/c. The differential 
spectra of k particles in pseudorapidity bins are 
analysed relative to the similar background spec- 
tra. The dips in the distributions of such ratios 



are found to be independent on the number k of 
particles, ~ 1.8. The number of clusters is es- 
timated to be about 2-3 with two particles per 
cluster in average. It is found that the data in the 
non-annihilation channel is similar to that from 
the inelastic pp collisions at 69 GeV/c from the 
same accelerator. These observations are treated 
in frames of different mechanisms, e.g. a model 
with mesons emitted from intermediate nucleon 
isobars is suggested. 

I am thankful to the ICHEP02 Organising Com- 
mittee, to convenors of the QCD Session, and to 
my colleagues at CERN and particularly in OPAL 
for giving me the opportunity to make this pre- 
sentation and for their kind help and support. 

REFERENCES 

1. I.M. Dremin, J.W. Gary, Phys. Rep. 349 

(2001) 301. 

2. V.A. Khoze, W. Ochs, Int. J. Mod. Phys. A12 
(1997) 2949. 

3. E.A. De Wolf, I.M. Dremin, W. Kittel, Phys. 
Rep. 270 (1996) 1. 

4. M. Siebel, Fragmentation of quarks and glu- 
ons, these Proceedings. 

5. S. Todorova-Nova, Bose-Einstein correlations 
at LEP and HERA, these Proceedings. 

6. P. Abreu, A. De Angehs, DELPHI 2002-052 
CONE 586 (2002), ICHEP02 abs. 232; 
DELPHI Collaboration, P. Abreu et al., Phys. 
Lett. B479 (2000) 118. 

7. OPAL Collaboration, G. Abbiendi et al., Eur. 
Phys. J. C19 (2001) 257. 

8. OPAL Collaboration, Physics Note PN511 

(2002) , ICHEP02 abs. 373. 

9. DELPHI Collaboration, P. Abreu et al., Nucl. 
Phys. B386 (1992) 471. 

10. OPAL Collaboration, G. Abbiendi et al., Eur. 
Phys. J. Cll (1999) 239. 

11. L3 Collaboration, L3 Note 2758 (2002), 
ICHEP02 abs. 494. 

12. Liu Feng, Liu Fuming, Liu Lianshou, Phys. 
Rev. D59 (1999) 114020 and refs. therein. 

13. OPAL Collaboration, G. Abbiendi et al., 
Phys. Lett. B523 (2001)35, ICHEP02 abs.319. 

14. E.G. Boos et al., ICHEP02 abs. 37. 



